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We present in analytic form the 0{a^) correction to the H ^ gg partial width of the standard- 
model Higgs boson with intermediate mass Mh < 2Mt. Its knowledge is useful because the 0{a'^) 
correction is sizeable (around 20%). For AIh = 120 GeV, the resulting QCD correction factor reads 

1 + (215/12)ai'^^ {Mh)/-k + 152.5 (af^ {Mh)/'A + 381.5 (af'> {Mh)/A ~ 1 + 0.65 + 0.20 + 0.02. 
The new four-loop correction increases the total Higgs-boson hadronic width by a small amount of 
order l%o and stabilizes significantly the residual scale dependence. 

PACS numbers: 14.80.Bn 12.38.-t 12.38.Bx 



INTRODUCTION 

Within the Standard Model (SM) the scalar Higgs bo- 
son is responsible for mechanism of the electroweak mass 
generation. It is the last fundamental particle in the 
SM which has not yet been directly observed. Its fu- 
ture (non-)discovery will be of primary importance for 
all the particle physics. The SM Higgs boson mass is 
constrained from below, Mh > 114 GeV, by experiments 
at LEP and SLC 1, 2]. Indirect constraints from preci- 
sion electroweak measurements J3| set an upper limit of 
200 GeV on Mh- 

Adopting the framework of the SM, the coupling of 
the Higgs boson to gluons is mediated by virtual massive 
quarks ^1 and it is this coupling which plays a crucial 
role in Higgs phenomenology. Indeed, with the Yukawa 
couplings of the Higgs boson to quarks being proportional 
to the respective quark masses, the ggH coupling of the 
SM is essentially generated by the top quark alone. The 
ggH coupling strength becomes independent of the top- 
quark mass Mt in the limit Mh ^ 2Mf 

The process of the gluon- fusion, gg — > H , provides a 
very important Higgs production mechanism over the all 
Mh range under consideration. The corresponding cross- 
section is significantly increased, by approximately 70%, 
by next-to-leading orderfNLO) QCD corrections, avail- 
able since long [a S S 13 • The largeness of the correc- 
tion along with the large residual scheme dependence of 
the result have motivated the calculations of the NNLO 
terms [3, llfl llH' Very recently even the leading N'^LO 
corrections to the inclusive cross-section have been com- 
puted fX2^ . As a result of these remarkable theoretical ad- 
vances the theoretical uncertainty of the production cross 
section is reduced significantly and is estimated around 



20%. 

The QCD corrections to the closely-related process — 
the production cross-section of the Hig gs d ecay into two 
gluons — are presently known to NNLO [ij, [iJI only. Re- 
cently it was pointed out in work |l5| that the ratio of 
the production cross-section to the decay rate is signif- 
icantly less (by a factor of 2) sensitive to higher order 
QCD corrections than the individual observables, since 
the corresponding K-factors are similar in size and tend 
to cancel to a significant extent. The work also argues 
that it is this ratio which presents the theoretical input 
to analyses of Higgs couplings extractions at the LHC. 
Thus, the knowledge of the N^LO QCD corrections to 
the Higgs decay rate into gluons is highly desirable. 

In this publication we present in analytic form the four- 
loop 0{al) correction to the H ^ gg partial width of the 
standard-model Higgs boson with mass Mh < "^Mt- 



CALCULATION AND RESULTS 

We start by constructing an effective Lagrangian, CcS, 
by integrating out the top quark |l3L llq . This La- 
grangian is a linear combination of certain dimension- 
four operators acting in QCD with five quark flavours, 
while all Mt dependence is contained in the coefficient 
functions. We then renormalize this Lagrangian and 
compute with its help the H ^ gg decay width through 
0(af). 

The effective Lagrangian can be written in the form 



-cff 



-2^/^G]^^HCi [0[] . 



(1) 



Here, [0[] is the renormalized counterpart of the bare 
operator O'l — G^'^^.G^''^" , where Gafj.u is the colour 



field strength. The superscript denotes bare fields, and 
primed objects refer to the five-flavour effective theory, 
Ci is the corresponding renormalized coefficient function, 
which carries all Mf dependence. 

Eq. ^ directly leads to a general expression for the 
H ^ gg decay width. 



TiH -. gg) = 
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where 



GG/ 2n 



n^"(<z 



e^'^^0\T{[O[]ix)[O[]m\0)d^ (3) 



is the vacuum polarization induced by the gluon operator 
at q^ = Mfj, with q being the external four-momentum. 
It is customary to write eq. (jJJ in the form 



r{H~.gg)^KrBorn{H^gg), 



(4) 



where {Gp is Fermi's constant) 

GpMfj 



rBorn(^ -^ gg) 
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and so-called iiT-factor reads: 
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The coefficient function Ci is known in N'^LO |l7lll8l 
l!?! and reads 
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Here 



-/j,t 



= 1-^ 



, with Mt being the on-shell top quark 
mass; we have displayed for generality the result with 
the effective number of light quark flavors denoted as ni . 
In the numerical evaluation we have set /z = Mi and 
ni = Uj — \ = h. Cn = C("-) is the Riemann's Zeta- 
function. 

Thus we are left with the calculation of the last factor 
in Eq. Q, namely the absorptive part 11'^'^ in N^LO, that 
is to a^. In fact it turned out to be more convenient to 
calculate the correlator © per se and take its absorptive 
part subsequently. Since 11'^'^ starts in the leading order 
from a one-loop diagram, the 0{al) calculation faces as 
many as 10240 four loop diagrams (at NNLO [13 the 
number was 403). 

The overall strategy of our calculations was identical to 
that used by us before in works |2(i,,2l| . First we generate 
the contributing diagrams with the package QGRAF |22| ■ 
Second, usin g th e criterion of irreducibility of Feynman 
integrals [23, |2J| , the set of irreducible integrals involved 



in the problem was constructed. Third, the coefficients 
multiplying these integrals were calculated as series in 
the 1/D — > expansion with the help of an auxiliary 
integral representation [23 ■ Fourth, the exact answer, 
i.e. a rational function of D, was reconstructed from 
this expansion. The major part of the calculations was 
performed on the SGI Altix 3700 computer (32 Itanium-2 
1.3 GHz processors) using the parallel version of FORM 
|26l I27I 1281] and took about two weeks of calendar time 
in total. 



After renormalization, our result reads 
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In order to better understand the structure of the af 
term in (|12(l it is instructive to separate the genuine four- 
loop contributions from the function 11'^'^ (q^) from essen- 
tially "kinematical" , so-called tt^ -terms originating from 
the analytic continuation. For a given order in as these 
extra contributions are completely predictable from the 
standard evolution equations applied to the "more lead- 
ing" terms in Il'~^'~^{q^) proportional to some smaller pow- 
ers of as . The corresponding expression for Im 11'^'^ as- 
sumes the form 



2q' 



imn 



GG 



1 + 12.4167 fls 
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(104.905- 36.257 ) g? 



037-1098.' 



where a'g stands for ai" /tt. 

In eqs. H9lllf) we have set /x^ = g^, the full /x depen- 
dence can be easily recovered with the standard RG tech- 
niques from the fact the product 



(/3(-)(a's))'n«««,M/g^ 



is scale independent (|la|)- In numerical form Imll'-^'^ 
reads (we set m = 5) 



2(7"* ^^ 

— i-Imn<^^ ^ 

TT 

1 + 12.4167 Os -I- 68.6482 a? - 212.447 a^ 



(12) 



where we have underlined the contributions coming from 
analytic continuation. Thus, the present calculation con- 
firms the pattern first observed on the example of the 
scalar correlator in work |i2l| : the "kinematical" tt^ terms 
tend to neatly cancel the genuine higher order contribu- 
tions. 

We are now in a position to find the 0{a^) term of 
the K factor in Eq. ^. To this end, we first multiply 
Ci by R'^{q'^ = Mjj), then eliminate al {jj) in favour of 

ai"'V) dlilE^I and, finally, choose the ^JL = Mh (to 
get a compact expression). The result reads: 



K = l + 17.9167a's -I- (156.81 
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If we also use Mt = 175 GeV, Mh = 120 GeV and 
of {M}j)/t^ = 0.0363 we arrive at 



K 



1 + 17.9167a's -I- 152.5(a's)2 -|- 381. 5(0^ 

1 + 0.65038 + 0.20095 -f 0.01825. (15) 



Thus, we observe that, unlike the NLO and NNLO cases, 
the 0{a^g) correction, being by a factor more than ten less 
than the previous a^ one, has quite a moderate size. The 
welcome stability of the perturbation theory is also con- 
firmed by testing the scale dependence of the i^-factor. 
By changing the scale (x from M^j^ to 2 Mu we find that 
the maximal deviation of the ii'-factor from its value at 
\x = Mh decreases from 24% (LO) to 22% (NLO), 10% 
(NNLO) and, finally, to only 3% at NNNLO! 



CONCLUSION 

We have computed the N^LO correction of order as to 
the H ^r gg partial width of the standard-model Higgs 
boson with intermediate mass Mh < "2 Alt- Our calcu- 
lation significantly reduces the theoretical uncertainty of 
the QCD prediction for this important process. 
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